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Abstract 



We study large deviations for some non-local parabolic type equations. We show that, under 
some assumptions on the non-local term, problems defined in a bounded domain converge with an 
exponential rate to the solution of the problem defined in the whole space. We compute this rate 
in different examples, with different kernels defining the non-local term, and it turns out that the 
estimate of convergence depends strongly on the decay at infinity of that kernel. 
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1 Introduction 

Consider continuous and bounded solutions u : M. N x [0, oo) — * R of the linear non-local 
equation 



where u(x, 0) = uq(x) is fixed, bounded and continuous. For simplicity, we will assume 
throughout the paper that solutions are non-negative, and thus also uq 0. The kernel J 
is assumed to be a symmetric, continuous probability density. 

The main contribution of this paper is to describe how solutions ur of Ql.ip , but defined in 
the ball Br = {x € K w : \x\ ^ R}, converge to the solution u of 

Let us first explain what ur is exactly: we consider here the notion of Dirichlct problem 
that consists in putting ur = not only on the topological boundary of Br, but also in all 
the complement of Br x [0, oo). In this way ur solves the equation 
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(1.1) 




(1.2) 




J (x - y)u R (y, t) dy - u R {x, t) in Br x (0, oo) 



with initial data ur(x,0) = uq(x) in Br. We refer to [1] and [5] for more information on 
these non-local Dirichlet problems, and also [2] for similar questions (with singular kernels). 

As one can imagine, under suitable assumptions, ur will reasonably converge to u as R — ► oo, 
but we want to obtain some estimates of how fast convergence occurs. Actually, this prob- 
lem may be seen as a numerical question since of course, computing numerically solutions 
requires a bounded domain. In this case one has to know how far from the real solution the 
computed one is. 

In the case of the Heat Equation, the answer is given in [T] : the distance between u and ur 
in the ball of radius 9R (with < 8 < 1) is estimated by 



which means that convergence occurs exponentially fast, with a rate of the order of R 2 
inside the exponential. 

Our aim is to produce similar estimates for non-local equations We face here several 

difficulties, which imply non trivial adaptations of ideas and techniques in pQ, that we list 
below: 

i) Various behaviours of J imply various rates: the importance of the tail of J enters into 
play, since the operator puts emphasis on the difference between u and ur = far from 
the point where we compute ur, as we shall explain below in the subsection devoted to 
the probabilistic aspects. Roughly speaking, the more J is big at infinity, the slower ur 
converges to u. 

ii) The structure of the Hamiltonian which describes the rate function X is completely 
different: for the Heat Equation the associated Hamilton-Jacobi equation is <9fZ+|VZ| 2 = 0, 
hence H(p) = p 2 . Here the problem for the rate function is related to the hamiltonian 



which, although it is local, is the limit as R — > oo of hamiltonians involving a non-local 
term, 



This localization process is one of the main interesting features of this problem. 

Hi) We are not facing here a diffusive effect, but more a transport effect: in the case of 
the Heat Equation, the scaling used in pQ in order to proof convergence is the parabolic 
(Rx, t) one (equivalent to (R 2 x, Rt)). Here, we have to use a hyperbolic change of variables 
(Rx, Rt) in order to scale the problem in a suitable way. 

Probabilistic context - The term "large deviation" comes from the french "grands 
ecarts" which was used first to describe how far from the normal distribution, some excep- 
tional events are. For instance, it is well-known that if (X n ) is a sequence of independent 
and identically distributed random variables with E(-Xi) = fx, then 




(1.3) 





X\ + X 2 + • • • + X n p 



ri 
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as n — > +00, the convergence occurring in law (this is the law of large numbers). Now, one 
may wonder how to estimate, for e > small, the quantity: 



Xi +X 2 



X n 



> e 



A result of Cramer (1938, see [7] for a proof), shows that if one defines the rate function 



1(e) 



lim 



log] 



X t +X 2 + 



x„ 



> e 



then 



which implies, 



/(e) = sup {et - logM(i)} , where M{t) = E[e tXl ] < oo . 
teR 



X\ + X 2 



x„ 



/' 



> e) 



-nl(e) 



This exponential behaviour is typical of what is called "large deviations" . 

In this paper, (u — ur)(-,T) mesures in some sense the total amount of process that can 
escape from the ball Br between times and T (see [T] and [I] for more explanations 
about this aspect). Thus, our results may be viewed as "large deviations" results in the 
sense that the probability of escaping the ball Br up to a given time becomes small as 
R — > oo. Exponentially small in fact, with a rate which depends on the tail of J since 
this tail measures the amount of "big jumps" . Values of J near the origin only concern 
"small jumps" that are not relevant as far as escaping the ball is at stake. So this is why, 
as we explain in Section 5, adding a singularity at the origin does not change the rate of 
convergence. 



Main results - After a preliminary section in which some properties of non-local and 
Hamilton-Jacobi equations are reviewed, Section [21 we devote Section [3] to the theoretical 
behaviour of the problem. Theorem l3.1l is a general result that gives the following estimate 



lim |u - u R \(x,t) < e -Ri~WWt)+oW 

R — >oo 

where the rate function Iqo satisfies the limit Hamilton-Jacobi problem associated to the 
Hamiltonian H defined in (11.31): 



| dJoo + ff(VIoo) =0 in Bx x (0, oo), 
(1.4) </oo = on dBi X (0,oo), 

I loo (x, 0) = +oo in B\. 

Using a Lax-Olcinik formula, see [8], we obtain a semi-explicit expression for /oo, sec (|3.1|) . 

In Section [4] we study different cases where explicit computations can be derived. A typical 
and interesting example, often considered by authors, concerns the case when J is compactly 
supported in, say, B v . We prove that the behaviour of loo is of "sins" type, which in turn 
implies that the rate function is of R\nR type; i.e. for any 6 € (0, 1), 



sup |u-Ufl|(t) <exp(-(l-0)^^ — WR + ^)r) as R 

\x\^6R V V ' 



OO. 



Other examples which imply different rates, like R(\nR)( a ~ 1 ^ a for J(y) ~ c~' y ' , are 
dealt with and the limit case J(y) = e~l y l is also considered even if this leads to a singular 
hamiltonian, defined only for \p\ < 1, a case not covered by the results of Section [3J 



3 



Finally, in Section 5 we explain how to extend these results to more singular situations, 
when the kernel is not integrable near the origin. Provided the decay at infinity remains 
reasonable (i.e., the Hamiltonian is defined everywhere), the presence of a singularity at the 
origin does not change the behaviour since as we have already mentioned, only the tail of J 
is important in estimating (u — ur) . 

Notice however that using kernels with singularities requires a suitable concept of solution; 
we use here the notion of viscosity solutions derived in [5] for Levy-type operators, which 
allows us to handle this situation. The Hamiltonian also needs to be modified a bit, as a 
corrector term appears in the equation: 

H(p)= [ (cPy-l-l {M<1} (p-y))j(y)dy. 

JR« v ' 

2 Preliminaries 

As a first step in order to prove the main theorem of this paper, Theorem 13. 1[ we need to 
state some properties of the non-local equations we are considering here, (jl.ip and (|1.2[) . 
We also need to do a thorough study, see Subsection 12.21 of the hamiltonian H and of the 
related Hamilton- Jacobi equation (|1.4[) . To this aim, let us first state the exact assumptions 
on kernel J: 

(2.1) J non-negative, radially symmetric, / J(y)dy = l. 

JR" 

The case J J(y) dy = c < oo can be also be considered just by doing a change of variables 
in t. 

These are natural conditions used to give sense to the non-local equation. But moreover, 
we assume that J has fast decay at infinity: 

(2.2) For any p e l w , H(p) = [ (c p v -l)j{y)dy<+oo. 

For instance, compactly supported kernels are authorized and non-compactly supported 
kernels also, provided they decay sufficiently fast at infinity. The limit case is J(y) = e _ ' y ' 
for which H is finite only inside the ball B\. We also give a formal derivation of estimates 
for this limit case, that we intend to prove in a forthcoming paper. Notice also that Section 
[5] extends the results to kernels J that are singular at the origin provided they remain Levy 
measures, and (|2.2j) is satisfied (but integrating on the domain {\y\ > 1}). 

2.1 Properties of the non-local equation 

Equations and (|1.2[) have been studied in [U [5| ; we refer to these papers for proofs 

of existence, uniqueness and comparison results, as well as other qualitative properties like 
positivity up to the boundary. 

However, we want to deal here with only bounded (and continuous) initial data, a case not 
covered in [1] since the argument of Fourier transforms that was used there requires that 
both uq and its Fourier transform are integrable. So let us first give a more general existence 
and uniqueness result in M. N . 

Definition 2.1 Let uo : R N — * K be continuous, non-negative and bounded. A strong 
solution of (jl.ip with initial data u(0,x) = Uq(x) is a function u £ C°(M. N x [0, oo)) such 
that Ut € C°(R Ar x (0, oo)) and ([l.ip holds in the classical sense, pointwise. 
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Before proving existence, we state a comparison result valid in the class of bounded strong 
solutions: 

Proposition 2.1 Let u,v be bounded strong solutions of with initial data uq and vq 

respectively. If uq ^ vq, then u ^ v everywhere. 

Proof. It is contained in [3J thm 3], using as Levy measure [i(z) = J(z)dz. The only 
adaptation is that we are here in a parabolic situation whereas the cited Theorem works for 
the elliptic case. The adaptation is standard, the "u t " term replacing the "— 7^" term in 
the Hamiltonian. ■ 

Theorem 2.2 Let uq : R N — > K be continuous non-negative and bounded. Then there exists 
a unique strong solution of (jl.ip with initial data Uq. 

Proof. Consider a sequence uq 71 S ^""(R^) n L 1 (M. N ), so that the Fourier transform of uq u 
is also integrable, and such that UQ n — > uq monotonically: one can choose a sequence of the 
form UQn = uq ■ Xn, Xn smooth, compactly supported and converging monotonically to 1. 
We know from [4] that there exists a unique solution u n of (11. ip which is in fact continuous 
since its Fourier transform remains integrable. 

As n increases, the sequence increases by comparison in the class of continuous and bounded 
solutions. Since (u„) is bounded, there exists a limit u defined in all M. N x (0, oo). The same 
happens with the sequence (v n ) defined as v n (x, t) := e t u n {x, t), which satisfies the equation 

d t V n — J * V n - 

Passage to the limit in the sense of distributions is done using dominated converge for the 
convolution term: J * v n converges to J * v and the limit equation dtv — J * v is satisfied in 
the weak sense. 

Now, using well-known properties of the convolution, we have that J * v is continuous (recall 
that J is integrable while v is bounded), so that d t v is a continuous function. The same 
holds for u and thus we recover a strong solution. 

Uniqueness follows from the above comparison principle, Proposition ^. II ■ 
Similar arguments show that indeed ur converges to u: 

Proposition 2.3 Let uo : R — > R be continuous non-negative and bounded. For any 
R > let un be the solution of (|1.2|) with initial datum Uq- Then the sequence (ur) 
converges monotonically as R — > oo to the solution u of (jl.ip with initial datum Uq. 

Proof. We use the same proof as that of Theorem 12.21 The sequence (ur) is increasing 
with respect to R and it is bounded by ||wo||- Hence it converges to some u which is a 
distributional solution of (jl.ip . Using again properties of the convolution, we get that u is 
in fact the unique strong solution of the equation with u{x, 0) = uq(x). ■ 

2.2 Hamiltonians and Lagrangians 

As it is well-known in the theory of Hamilton- Jacobi equations, the Legendre-Fenchel trans- 
form L of H defined as 

(2.3) L(q) = sup {p ■ q - H(p) } , qeR N . 

pGR N 
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plays an important role in representing solutions. Since the rate function 1^ is the solution 
of a Hamiltin- Jacobi equation, see (|1.4p we review first some basic properties of both H and 
L. 

Lemma 2.4 Let J satisfy (|2.1[) and (|2.2p . Then both H and L are nonnegative, increasing, 
strictly convex, superlinear and radially symetric. 

Proof. In order to prove rotation-invariancy, let a be a rotation of the sphere S N_1 , then 

H(ap) = [ (c^>y~l)j(y)dy 

( e P(-- 1 y)-l)J(y)dy 
(e™' -l)j(c-y')dy' 



c 



p-y 



l)J(y')dy' 



Hip) 



(remember that J is symmetric). 

Strict convexity and superlincarity come from the estimates: 

D 2 H(p) = f cP-y\y\ 2 J(y) dy > , D p H(0) = / yj(y) = . 



Indeed, by rotation-invariancy, D p H(p) is always pointing in the direction of p 
and moreover 

(2.4) p.D p H(p) = \p\\D p H(p)\ >0 

for any p ^ 0. Using the strict convexity of H, we get that H grows faster than linearly 
along the lines Xp, A > 0. 

Strict convexity, together with H(0) = imply that H is nonnegative. 

Finally, it is well-known, [5], that if H enjoys the above properties, so does L. ■ 

In the case we are considering, i.e. J radially symmetric, since H and L are also radially 
symmetric, throughout the paper we denote by H* and L* the functions defined on R+ such 
that 

(2.5) H*(\p\)=H(p) and L* (\q\) = L(q). 

The following technical trick will be used several times in the paper, hence we state it as a 
Lemma: 

Lemma 2.5 Let J satisfy (|2.ip and (|2.2I) . For q fixed, let po = po{q) be a maximum point 
defining L(q), that is, L(q) = p ■ q - H(p ). Then p Q ■ q = \p \\q\. 

Proof. Notice first that by definition of L, po solves the equation q = D p H(po). The result 
follows using (|2.4p . ■ 



G 



2.3 Viscosity solutions for the rate function equation 

In Section O wc will have to study the following Hamilton- Jacobi equation with Cauchy- 
Dirichlct boundary values, 

!d t I A + H{VI A ) = inSx x (0,oo), 
I A = on dBx x (0, 00), 

I A {x,0)=A in Si. 

Let us first recall the definition of viscosity solutions for this equation (see for instance [6]): 

Definition 2.2 A locally bounded u.s.c function u : B\ x [0, 00) — > R is a viscosity subso- 
lution of (|2.6j) if for any C 2 -smooth function tp, and any point (xo,to) € B\ X [0, 00) where 
u — if reaches a maximum, there holds, 

xqEBi => d t ip + H(Vp) < at (x ,t ), 
xo € dB\ =>• min {dtf + H(Vtp) ; u} ^ at (xq, to) , 
t = u^A in B\. 

A locally bounded l.s.c. function is a viscosity supersolution if the same holds with reversed 
inequalities and min replaced by max at the boundary. Finally a viscosity solution is a 
locally bounded function u such that its u.s.c. and l.s.c. cnvcloppcs arc respectively sub- 
and super-solutions of (|2.6[) . 

Since H is convex we have the following representation: 

Lemma 2.6 If J satisfies (|2.1|) and (|2.2p . the unique viscosity solution of ()2.6p is given 

by, 

I A (x,t)=min{A,tL*( dist{x ' dBl) )} m B lX (0,oo). 

Proof. Recall that the assumptions on J imply that both H and L are finite everywhere, 
convex, radially symmetric and super-linear, see Lemma [2. 41 We then start from the Lax- 
Olcinik formula in the bounded domain B\ x (0, 00), see [5], 

I A (x,t) = min (tL(y^) + A] A min ! ( t - s)l(^-)X , 

(we denote by a A 6 the infimum of a and b). Since L is symmetric, we can rewrite it using 
L* . The fact that L is nonnegative and L(0) = implies: 

(2.7) I A {x,t) = AA min Ut - s)L* ( fcM) j . 

Since L* is increasing, the min is attained at the point y such that \x — y\ = dist(x, dB\). 
Notice now that since L is convex and L(Q) = 0, then for any fixed c > 0, and variable 
r > 0, 

{L *y {cr)> EM. 

cr 

This implies that the function r i— > L*(cr)/r is increasing so that the minimum in (|2.7[) is 
attained for s = (use c = dist(x, dB\) and r = (t — s) _1 which is minimum for s = 0). 
Combining all these estimates, we are led to Lemma \2. 61 ■ 
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3 Theoretical Behaviour 



The main goal of this section is to derive a theoretical bound, in terms of the Lagrangian 
L, for the error made when approximating the solution, u, of by solutions, Ur, of the 
Dirichlct problem (11.21) . 

Theorem 3.1 If J satisfies (|2.1[) and (|2 .2[) . then for any fixed x € and t > i/iere 
holds 

Urn |«-u fl |(ar,*) < e -*M*/*.t/*)+o(i)H 

i? — >00 

where the rate function is given by 

(3.1) !„(*,*) =tL*(**^2). 

Moreover, for any 9 € (0, 1), < ti < t2 < +oo i/ie /zmii is uniform in the set 

{\x/R\ <0, t/R€ (ti.ta)}. 

3.1 The transformed equation 

Let us denote by ^ = « - the solution of (|1.2p in S^ x (0, oo) with initial value 
vr(x,0) = for x £ Br and "boundary data" Ur = u for |a;| 5^ R. 

Now we first rescale the equation both in x and t as follows: 

wr(x, t) = vr(Rx, Rt) for x € Si, t . 

Then wjj satisfies a rescaled equation in the fixed ball B\, with rescaled nucleus Jr(x) = 
RJ(Rx): 

dtW R (x,t) = R[(Jr * w R )(x,t) - w R (x,t)] . 
Let K = K(uo) = ||uo||oo and consider ip the solution of the following Dirichlct problem: 

d t ipR.(x) = R[(Jr * i/>r)(x) - iPr{x)] in B x x (0, oo), 
< ip R {x, t) = K in (R N \ B{) x (0, oo), 

>i?(x,0)=0 in Si. 

Since ^ Vr ^ K , then a standard comparison yields wr ^ ip R . 

In order to estimate tpR we follow [1] and perform the "usual" logarithmic transform, but 
we have to rescale accordingly, dividing by R (and not R 2 as it is the case for the heat 
equation). So, remembering that for t > 0, ipR > 0, let us define 

Ir(x, t) = -— ]xl(i/} R (x, t)) . 

Then 

d t ^ R {x,t) = -Re- RI *^d t I a (x,t) 

and 

jR*ipR(x,t)-ijjR{x,t) = / J R (x-y)(il;R(y,t)-ipR(x,t))dy 

JR N 



[ RJ(Rx - Ry)e- RlR{x ' t) ( e fl i z *<»'*W*0M» _ i) dy, 
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which becomes, if we do the change of variables y — x — s/ R, 

c -RI n (x,t) f j/ a \rR{I R (x,t)-I R (x-s/R,t)}_- i \ ds ^ 

We arrive at the following equation for Ir: 

d t I R (x,t) = f L-n{^+v/R,t)-^,t)} A J{y)dy: 

JR N V ' 

which formally converges to the Hamilton- Jacobi equation, 

(3.2) d t I + H(X7I) = with Hip) = J {c p - y - l)J(y) dy . 

To justify convergence of Ir towards the solution of (|3.2p , we have to use viscosity solutions. 
This is done in the next subsection. 

3.2 Limit Hamilton-Jacobi Equation 

Thanks to modulus of continuity estimates proven in [5], and the fact that u, ur are bounded, 
we could extract a subsequence that converges locally uniformly; but we shall however 
use the "half-relaxed limits" method to handle the hamiltonian. 

A first problem comes from the fact that if ipR approaches zero, then Ir may not remain 
bounded. Hence to avoid upper estimates for Ir, we use the same trick as in [T] which 
consists in modifying Ir a little bit. For any A > 0, let 

l£{x,t) = ~\n{i> R {x,t)+e- RA ), 

which is bounded from above by A. Let us notice that since equation (jl.ip is invariant 
under addition of constants, Ir satisfies the same equation as Ir. 

Proposition 3.2 The sequence {Ir) converges locally uniformly in B\ x (0,oo) as R — * oo 
towards the unique viscosity solution I A of (|2.6p . 

Proof. We introduce the half-relaxed limits, 

I A (x,t) := limsup*/^(a;,i) = limsup Ir(x' \t') 

R^oc- («',t')->(*,t) 

-R-.00 

and 

L A (x,t) := liminf *Ir{x,£) = liminf l£(x',t'), 

R~ +00 (x',t')^(x,t) 

H-.OCJ 

and we shall prove that they are respectively viscosity sub- and super-solutions of the limit 
problem (|2.6p . Then a uniqueness result will allow us to conclude. 

^4 

Let us take a test function ip such that / —tp has a maximum at {xq, to). Up to a standard 
modification of ip, we can assume the maximum is strict so that there exist sequences 
R n — > +oo and {x n ,t n ) — > (xo,to) such that 

l£ — ip has a strict maximum at (x n , t n ) . 
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Case 1: the point (xo,to) is inside Bi x (0, oo). Then for n big enough, all the points 
(x n , t n ) are also inside B\ x (0, oo) so we may use the equation for Jj£ at those points and 
pass to the limit. 

Since dtl^ n is continuous, we have dtl^ = dt<p at (x n ,t n ) and moreover, for any z £ R*, 

^R n i x n z >tn) ~ ^R n {Xni tn) ^ fiXn + Z, t„) — (p(x ril t n ) . 

Using this, we fix e > and split the equation for l£ into two terms as follows: 



dMx n ,t n ) < -/ ( e -nW(^v/Rn,t n )-^ n ,t n )}_ 1 \j (y)dy 

J\y\<M V ' 

+ / {lpR n (Xn +y/Rn,tn) ~ 1pR n { x n,t n )}j(y)6 

J\V\>M 



Since ipn is bounded by K(uq), wc can choose M big enough so that the second term is less 
than s, independently of n. 

For the first term, we write a Taylor expansion for ip near point x n : there exists a € Bm 
(the ball of radius M) such that 



£ . 



y\<M 



Since remains in Bm and ip is smooth we have that _D 2 </?(£„) remains bounded. Hence, 
we can pass to the limit asm +oo: 



d, 



<p(x , to) < - / (c-WM) -v - l) J(y) dy + e 

J\v\<M V ' 



\y\<M 

Since H(\7(p(xo, to)) < +oo, we can let e — > and M — > +oo to get in the limit 
d t <p{x , t ) + f f c -^M)-y _ A J{y) dy <; o . 

JR N K ' 

A 

This shows that I at (xo,to) is a subsolution. 

Similar calculations lead to the supersolution condition at (xo,to) for I . 

Case 2: the point (x ,to) is located at the boundary, xo e dB\. Then the sequence 
(x n , t n ) may either lie inside B\ x (0, oo), or we may have \x n \ 1. In this last case is where 
the relaxed boundary condition in the viscosity sense, sec Definition 12.21 comes from. 

If x n € B\, we may use the equation as in the previous case, and we can do so even if 
x n G dB\ since for Jj| the equation holds at the boundary (see [5]). If on the contrary 
\x n \ > 1, then In n {x n , t n ) = so that in any case, one has 

min {d t p + H(Vip) ; l£ n } s$ at (x n , t n ) , 

and we pass to the limit as n — ► +oo to get the relaxed condition for / at the boundary. 
The converse condition for is obtained by the same method, with reversed inequalities. 

Conclusion: Using comparison between usc/lsc sub/super solutions for (|2.6p . we get 

the inequality I_ A ^ I , which implies equality of both functions. Hence, all the sequence 
converges to the unique solution I A . ■ 
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3.3 Proof of Theorem l3~7D 

This result only comes from the fact that for any A > 0, by construction 

1 A = inf (I, A), I A = inf (I, A), 



with 



I{x,i) := limsup *Ir(x, i), L{ x >t) := liminf *In(x,t). 



The fact that / = / , together with Proposition 12.61 yields the result, passing to the limit 
as A — + co. 



4 Estimates of the Lagrangian - explicit bounds 

Our next aim is finding estimates for the behaviour of L(q) as \q\ — > +oo, which in terms 
means estimates for 1^ . Although we consider here different J's the sketch of the proofs is 
always the same: 

1. At the maximum point, poj of L(q) it holds D p H{po) = q. 

2. Find estimates for D p H(p). 

3. Show that poq ^> D p H(po) and hence L(q) ~ p$q. 

4.1 Compactly supported kernels 

Let us begin with an explicit 1-D example: 

J(x) = 2 1 {kl<i}( x ) ■ 
For this choice of J, a straightforward calculus gives H(p) = (sinhp)/p, so that 

r sinhp-. 
L(q) =sup {pq j. 

peR P 

At the maximum point po = Po(q), we have 

coshp sinhp 

(4.3) q = — , 

Po Po 

gPO 

and since g^oowe necessarily have po — > oo. Hence (14.31) is equivalent to as q — > co. 

This implies, taking logarithms, that 



lng = p (<?)(l - f ) ^Po(q) 



so that as q — > +co, we have, using that sinhpo/po ~ 9, 

(4.4) ~glng-^^ + l~glng. 
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Remember that L is symmetric so that the same behaviour holds as q — ► — oo. 

This result can be easily extended to several dimensions just by using the symmetry of H(p). 
Without loss of generality we may assume p = Aei, which leads to a 1-D problem. 

In fact this "glng" behaviour is representative of what happens in general for compactly 
supported kernels in R N : 

Lemma 4.1 Let J satisfy (|2.1[) , (|2.2[) and supp ( J) = B n . Then the following behaviour 
holds: 

t*(\ w M ln M , , 

L ml) as \Q 

V 

Proof. Let po — Po{q) be the point where pq — H(p) attains its maximum. In order to 
estimate L(q) we have to investigate the behaviour of q — D p H(p) at Pq\ 

p-D p H(p) = [ p-yePyj(y)dy^ f \p\ V e^J(y)dy = \p\r ) e^. 
Taking logarithms and dividing by rj\p\ we conclude that 

it r\ r Hp-d p h( p )) f Hv\ P \) v\ P \\ 

(4.5) hmsup 1— | < hmsup j— ; 1 — |— , = 1. 

ipHoo vm ipi->oo V mp\ mv\) 

In order to obtain an lower bound for p ■ D p H(p) we split it into two integrals as follows: 
p-D p H(p)= f p-yePyj(y)dy+ f p ■ yf* J(y) dy. 

J{py^0} J{p-y>0} 

The first one is bounded from below by — \p\rj, using that J B J(y)dy = 1. For the second 
one we define for e > and a < 1, the set 

eta = {y' 1 <\y\<v,P-y>0-- e)\p\\v\ > «}• 

a 

Since H is radially symmetric and increasing with then 

D p H{p) = \{p)p for some \{p) > 0, 
and p ■ D p H(p) has its main contribution in Cf a . Hence 

p ■ yeP yj(y) dy > \p\\y\(l - e) [ e^^J(y) dy 
{p-y>0} Jc+ a 

j> |p|^(l_ e ) e a-=)H3 / J( y )dy 



^ o^e, a)\p' 71 

Summing up 



> C(e,a)\p\-(1- e)c^-^ p ^ . 
a 



p ■ D p H{p) > C( £ ,q)|pK(1 - £ )c (1 - £)|p| ° - \ p \t, > KC(e,a)\p\±{l - e)c (1 - e)|p| - , 
a a 

for some constant K. Therefore, arguing as in (14. 5[) we obtain for every a and e 



hminf HP ' DpH{p)) > hminf ( lnC ^ a ^ , MSN(l-e)) , *IpI(1 " ^ _ 1 " e 



IpI-oo 7/|p| IpHoo V »7b| ? /bl / " 
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Now, letting e —> and a — > 1 we conclude 

(4.6) MP^Ml „ !, 

v\p\ 

as |p| — » oo. The main point now is that D p H(p) = \D p H(p)\-£r which implies that (|4.6|) 
becomes 

1 _ ln|£> p g(p)|+ln(|p|) \n\D p H(p)\ 

v\p\ v\p\ 

At the maximum point q = D p (H (po)) we then have: 

\q\ln\q\ 



Pa ■ q = \Po\\q\ 



V 



The final step consists in proving that H(po) is negligible compared with po ■ q. To this aim, 
let Ji = II {lyKjji}) for i = 1, 2 so that J\ < J < J<x (remark that j Ji ^ 1). Thus, from the 
computation done before, we obtain 

These inequalities together with (|4.4[) yield the desired result. ■ 

We are now ready to give our first concrete estimate for compactly supported kernels: 

Corollary 4.2 Let J satisfy (|2.ip . (|2.2[) supp ( J) = £>,,. Then the following estimate 
holds: 

( (1 - /(I - 0)R\ , , \ 
sup lu-UflK^^exp — In — + o(l)R) as R -> oo. 

|a:|<0ii V ?7 V t / / 

Proof. According to the scaling (i?x, Rt) described in of the Proof of Theorem l3.11 we have 
to estimate 

j / x t \ t f&ist(x/R, dB{)\ t /dist(x, DBr) 
^KR'RJ^R V t/R ) ~ R V t 

To this aim let us keep t > fixed and take x € Res. for some 6> S (0, 1), so that 
dist(x, dBf>) ^ (1 — 6)R — > oo. We use Lemma |4~T1 to get that, as R — > oo, 



\R : RJ R \ t J 



(1-60-R, /(l - 6>)J2 
In ' 



77 V t 

We thus obtain the bound in BgR\ 

sup |„ - tta |(t) < exp ( - <I^>* In (<±ZM) + o{1)R ) , 

or in a simpler form, if t remains bounded: 

sup \u - u R \(t) < exp C - — — — In i? + o(i? In R) 

\x\^8R ^ V 
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We thus have a convergence rate of order R\nR. ■ 

This shows that we are not in the case of the heat equation for which the order is R 2 ; here 
the convergence occurs at a slower speed which is due to the fact that more paths of the 
process escape from the ball Br. But nevertheless, we get a somewhat good control of the 
error. 



Hence 



4.2 Fast exponential decay 

Lemma 4.3 Let J(y) = e - '^' , for a > 1. Then the following behaviour holds: 

as |«|-+oo. 

Proof. We shall do the calculations in 1-D, the adaptation to several dimensions follows 
the same lines as for the case of compactly supported kernels. 

We have 

H(p) = [ (e p - y - l)J(y)dy = [ e p y -^" dy. 
Jr Jr 

D p H(p) = f yffV-M" dy 

ye w-\v\ a dy + [ ye p - y -\ y \ a dy. 
{py>o} J{py<0} 

The integral over {p ■ y < 0} is bounded by some constant independent of p since in this set, 
e p ' y is less than 1, hence we will neglect it in the following estimates. 

For the integral over {p ■ y 0} consider the case p > 0, p — > +oo (the case p — > — oo being 
similar), and hence y > 0. let yo(p) = (p/a) 1 ^ a ~ 1 \ the point where p ■ y — y a attains its 
maximum. Then, since p ■ y — y a is non-decreasing in (yo — 1, yo) and a > 1 we get 
r rVo 

/ ye p - y - ya dy^ ye p - y - y " dy ^ (y - l) e p - ( - y °~ 1) ~ ( - Vo ~ 1)a > (y - l^™'^ e~ p . 

Taking logarithms, 

\n(D p H(p)) > In ( - 1) py - _ p 



pyo - Vo pyo - 2/o PVo ~ Va pyo - 2/o 

and hence replacing 2/0 by (p/a) 1 ^" -1 ^ 

,. . „ \n(D v H(p)) 
liminf - \ p y^' ' > 1, 

where 

On the other hand, let y\ = (p/Aa) 1 ^" ^ the point where p-y — \y a attains its maximum. 
Then, for A < 1 

ye p y-y a dy = ( ye p y - Xya e {x - 1)ya dy 
= C X G pyx - Xy ". 
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Since p ■ y\ — Xy" = c\(a)p a ^ a 1 \ with c\{a) = X 1 '^ a ^c(a), we have taking again 
logarithms and dividing by c(a)p Q / ( - Q-1 ), 

]n{D p H(p)) < c x (a) , ln(C A ) 



c(a)p«/( Q; - 1 ) ^ c(a) c(a)p Q /( Q - 1 ) ' 
and hence, for any A < 1, 

,. HD P H(p)) . c x (a) 
urn sup ~ \ / 1 TT ^ / \ • 

p-»oo c(a)p Q /( Q - 1 ) c(a) 

Letting now A — ► 1, we obtain 

hmsup , , ty ^ 1, 

and thus it turns out that 

ln(\q\)=\n(D p H{p))^c(a)p^ a - 1 l 

This estimate for D p H(p) yields that L(q) behaves as |(?|(ln |g|)^ Q_1 ^ Q provided we show 
that H (p) is negligible. 

To achieve this last step, let us prove that H(p) is at most of the order of D p H(p). Wc 
separate the integrals in two terms as follows: 



e p-y-\y\ a dy 

} 



H{p) = ( e p-y-\yr d y + [ 

J{\y\<i} J{\v\>i 

2 p + f ye p y - M " dy ^e p + D p H(p) 



< e p 

2D p H(p), 

(of course, this is valid for \p\ large). Thus the lemma is proved since 
L(q) = p Q q - H(p ) ~ p Q q ~ |g|(ln \q\)^l a . 



4.3 Critical exponential decay 

In this section we consider a case not covered by the results of Section [3] since we study the 
case J(y) = e~' 1 '', which is critical. Indeed, the corresponding Hamiltonian is only finite for 
H < 1: 

H(p) = U e ™-Mdy={^ forW<1 ' 
2 J r n oo for \p\ ^ 1. 



The supremum in 



L(q) = sup {p ■ q - H(p)} = sup {pq - 2 - j 

pen \ P \<i 1 -V > 



4p 

is obtained for q = =-r. Here, q — > oo corresponds to p — > ±1 which means that 

(1 — p 1 ) 2 

\l~p 2 \-2\q\"\ 
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so that 

However, Theorem 13.11 does not apply as such here, since H takes infinite values. So we 
only give a formal estimate for the sake of comparison: 

\u-u R \(x,t)^e-^ R+ °^ R . 

We shall adress the details of this case in a forthcoming paper, which imply non trivial 
adaptations of Section [3l 



5 Generalization to infinite activity jump diffusions 

In this section, we briefly explain how to extend our results to a class of singular kernels. 
We consider here functions J satisfying: 

/ fl A \y\ 2 ") J(y) dy < +00 , 

that is, we are interested in Levy measures with density J. Typical examples of such Levy 
measure are: 

1 e~ |y| 1 

where < a < 2. The first example is related to the well-known fractional Laplacian, while 
the second example is called "tempered a-stable law" among the probability community. 
The third example is singular at the origin, but compactly supported. 

However, the equation has to be understood in a special way: since J is not integrable near 
the origin, the equation should contain an extra term (called "corrector" ) in order to give 
sense to the integral term: 

d t u= / lu(x + y) — u(x) — 1 {|j/|<i}(Vii • y)\j{y) dy . 

Of course, if u were C 2 -smooth, the integrand would be close to D 2 u for y small, and 
everything would be integrable. But since we do not know a priori the regularity of u, we 
have to replace it by some smooth test-function and this is where viscosity solutions enter 
into play. We refer to [5J [5] for precise definitions and properties of visosity solutions in 
presence of Levy-type non-local terms. 

With this tool, everthing works exactly as we did for the non-singular case, except that the 
new Hamiltonian also involves a corrector term: 

(5.7) H(p)= [ {c?y -l-(p-y)l {lyl<1} }j(y)dy, pel". 

Notice that integrability near the origin comes from the assumption on J (since c v y — \—py ~ 
(py) 2 /2 for y ~ 0), but we have to face the same integrability condition at infinity: we shall 
assume that H(p) < +00 everywhere in M. N . 

So, the case of fractional Laplace operators is not covered here since in this case, H = +00; 
but the third example above can be dealt with. In view of Section l4~3l we hope also to cover 
the a-stablc law case soon. 
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For these adaptations, using the fact that the various manipulations that are used in Section 
[3] are valid for viscosity solutions, it can be checked that the same limit Hamilton- Jacobi 
equation is obtained: 

d t I oo +H{VI oo ) = 0. 

Now, when one takes a look at the singular Hamiltonian H, it appears clearly that the 
behaviour at infinity (i.e., for \p\ large) does not depend on the corrector term which is of 
lower order. Another way of understanding this is to think in terms of the total amount of 
process that can escape the ball Br: small jumps are not responsible for that behaviour (in 
first approximation). The main important contribution comes from big jumps, related to 
the tail of J. 

To illustrate this heuristic remark, let us consider the case of a compactly supported measure 
with a singularity: we consider the 1-D kernel 

J{y) = |3T+^ I [-i,i](y)> where «e(o,2). 

We do the same computations with 



which gives on the one hand D p H(p) ^ c(a)e p , so that 

lim sup — r-?- — ^ 1 , 

p^+oo c(a)eP 

and thus the maximal behaviour of L(q) is: 

hmsup — — — < c(a). 
p^oo \q\ln\q\ 



On the other hand we use the fact that the Lcgcndre-Fenchel is non-increasing, which is 
obvious from the sup formula relating L and H : since 

/ 1 e py ^ r 1 e py 

we know that the corresponding L satisfies: L(q) ^ L(q). But the modified kernel is a 
smooth and compactly supported function for which the behaviour is known: 

L(?)~</(a)|9|ln| g |. 

Finally we get that for some constants c\,C2 > 0, 

ci|g| In |q| < L(q) < c 2 |g|ln|g| . 



The conclusion is that the presence of singularities at the origin does not change the scale 
essentially Similar examples can be derived for non compactly supported kernels. 
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6 Numerical experiments 



We gather now some numerical examples, with the aim of illustrating the convergence 
theorem 13.11 

We try to approximate solutions of (|1.1[) . with initial datum uq = 1. The advantage of 
taking this datum is that the solution of the problem is u = 1. We have an exact solution 
that we do not compute numerically. 

In order to approximate ur we use a fixed mesh scheme that discrctizes the interval [-R, R], 
and an ODE integrator provide by Matlab® to integrate in time up to the fixed time 
to = 0.1. For each time-step, the integral involved in the term J * ur is approximated by 
the classical trapezoidal rule. 

6.1 Compactly supported kernels 

We consider here a compactly supported nucleus 

J(y) = ~(x-2)(x + 2). 
In Figure Q] we show convergence of ur to u — 1 . We observe that ur increases as R 




Figure 1: J compactly supported 
increases (the values of R are R = 10, 15, 20). 

6.2 Gaussian example 

Consider the case 

J(y) = -^=c- y2/2 . 
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Using the same ideas as in Section 0] we compute 

H(p) = -^L / c pv - y2/2 dy- 1 



2tt 
1 



'2tt 



-(p-y) 2 /2. eP 2 /2 d2; _ 1 



Now, D p H(p) = pe p I" 2 and is obtained for q = pc p I 2 . For q — > oo, we have ln< 
lnp + p 2 /2 ~p 2 /2, thusp~ 2(lnq) 1 / 2 so that 

L( (? )^2 g (ln< Z ) 1 /2. 



In the picture on the left of Figure [5] wc show again convergence for different values of R 
(R = 10, 15,20). The picture on the right represents for x and t fixed the convergence of 
this point to the exact solution. 



1.0004 




Figure 2: J gaussian, R = 10, 15, 20 



6.3 Critical exponential decay 

Finally we consider the case of a nucleus with critical exponential decay, 

As we mentioned before, this case is not covered by the results of Section [31 but it shall 
illustrate how convergence works for a general case: 



The picture on the right shows convergence of a point that moves with R; i.e. we take a 
point of the form (OR, t) with t = 0.1 fixed and 6 = 0.8. 
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12 14 16 18 20 



Figure 3: J critical exponential decay, R = 10, 15, 20 
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